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Abstract 

We derive the pairwise peculiar velocity distribution function of dark matter par- 
ticles applying the dark matter halo approach. Unlike the previous work, we do not 
assume a Gaussian velocity distribution function of dark matter in a single halo, but 
compute it self-consistently with the assumed density profile for dark matter halo. 
The resulting distribution function is well approximated by an exponential distribu- 
tion which is consistent with the previous observational, numerical and theoretical 
results. We also compute the pairwise peculiar velocity dispersion for different den- 
sity profiles, and provide a practical fitting formula. We apply an empirical biasing 
scheme into our model and present prediction for pairwise peculiar velocity dispersion 
of galaxies, and reproduce the previous results of simulations using our semi-analytical 
method. 
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1. Introduction 

Since Davis & Peebles (1983) first analyzed the anisotropy in the galaxy distribution 
from Center for Astrophysics (CfA) redshift catalog, it has been recognized that the peculiar 
velocity field of galaxies induces a significant systematic effect in statistics of observed galaxy 
distributions in redshift space. In particular, virialized random motion of galaxies produces 
an elongated pattern of galaxy distribution along the line of sight, called finger- of- God. This 
effect significantly suppresses the amplitude of the two-point correlation function of galaxies in 
redshift space, especially on scales below ~ 3/i _1 Mpc. 

The proper account of this redshift-space distortion requires a detailed model for the 
pairwise velocity distribution function (hereafter, PVDF) of galaxies. Davis & Peebles (1983) 
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discovered that the PVDF of the CfA galaxy sample is approximately described by an expo- 
nential distribution, instead of a Gaussian: 



where r±2 is the separation length and vu denote the pairwise peculiar velocity between a 
pair of galaxy along the line-of-sight direction. The quantity 0-12(712) is the peculiar velocity 
dispersion (PVD). The exponential form of PVDF was confirmed also later by analyses of N- 
body simulations of dark matter particles and of other samples of galaxies (e.g., Efstathiou et 
al. 1988, Fisher et al. 1994, Marzke et al. 1995). 

Theoretical models for the origin of the exponential PVDF of dark matter were put 
forward by Sheth (1996) and Diaferio & Geller (1996), and more recently by Sheth & Diaferio 
(2001). They phenomeno logically introduced a nonlinear model of PVDF using the Press- 
Schechter formalism. To be more specific, they assume that any dark matter particle belongs 
to one of virialized clumps (dark halos) with the 1-point velocity distribution function being a 
Maxwellian form. If one considers sufficiently small scales, the particle pairs of those separations 
are likely to be in the same halo, and then their PVDF is approximately given by 



where n(M) is the mass function of dark halos, /i2,ih(^i2|-^0 is the PVDF of dark matter 
particles within a halo of mass M. The quantity N p3i i T (ri2\M) represents the statistical weight 
proportional to the number of particle pairs with separation r i2 in the halo: 



where p(r\M) is the density profile of the halo of mass M, and So is the Dirac delta function. 
Adopting the singular isothermal distribution as a particular choice of the dark halo profile, 
Sheth (1996) showed that the scale-free model of P(k) oc k n with n = — 1 exactly reproduces 
the exponential PVDF (fj). 

While a perturbation theory (Seto & Yokoyama 1998; Juszkiewicz et al. 1998) also 
qualitatively explained why the Gaussian initial models approach the exponential PVDF, the 
above model is much more successful quantitatively. Further, a significant influence of the 
finger-of-Got effect appears at small scale, where the perturbative approach cannot be applied. 
Therefore in the present paper, we attempt to improve the Sheth (1996) model (|2|) for the PVDF 
in several aspects; first, we consider more popular CDM models instead of the scale- free power- 
spectra. Second, we adopt a series of more realistic density profiles for dark halos (Hernquist 
1990; Navarro, Frenk & White 1997; Fukushige & Makino 1997, 2001a; Moore et al. 1998; Jing 
& Suto 2000). Third, we derive the one-point PVDF of dark matter particles in a halo directly 
from the Abel integral of the above density profiles, instead of assuming the Maxwellian form 
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a priori. This approach is important since we can incorporate the scale- and mass- dependence 
of the PVD in a consistent fashion unlike the previous modeling. Finally we also apply the 
selection function following Jing, Borner & Suto (2002) so as to phenomenologically attempt 
to predict the PVD of galaxies out of that of dark matter particles. 

This paper is organized as follows; Section 2 describes our improved modeling for the 
PVDF on the basis of the dark matter halo approach. In §3, we present the resultant PVD in 
various cosmological models and discuss how the underlying halo profiles are sensitive to those 
results. We also provide a simple fitting formula of the PVD in the currently popular spatially- 
flat CDM model, which is useful in modeling the redshift-distortion effect. Then we attempt to 
consider the effect of the spatial biasing of galaxies relative to the dark matter particles on the 
PVD by applying a phenomenological biasing scheme. Finally section 4 is devoted to summary 
and conclusions. 

2. A Dark Matter Halo Approach to Compute the Pairwise Velocity Distribution 



will describe the details of the procedure below. Throughout the paper, we consider the three 
representative CDM models parameterized by the density parameter Qq, the dimensionless 
cosmological parameter X , the amplitude of the mass fluctuation smoothed over the top-hat 
radius of 8/i _1 Mpc, as, and the Hubble constant in units of 100 km/s/Mpc h; standard CDM 
(ft = 1.0, A = 0,<r 8 = 0.6,fc = 0.5; SCDM), lambda CDM (O = 0.3, A = 0.7, a 8 = 1.0, h = 0.7; 
LCDM), and open CDM (ft„ = 0.45, A = 0,a 8 = 0.83, h = 0.7; OCDM). Those models are 
normalized to satisfy the X-ray cluster abundances (Kitayama & Suto 1997). 

2.1. Density profile 

The density profile of dark matter halos plays a key role in our method. Following the 
recent suggestions from high-resolution numerical simulations, we adopt the following specific 
form: 



In the above, M is the mass of the halo, p(z) = f2op c o(l + z) 3 is the mean density of the universe 
at z, p c0 is the present critical density, S C (M) is the characteristic density excess, and r vir (M) 
and r s (M) indicate the virial radius and the scale radius of the halo, respectively. 
The virial radius is defined according to the spherical collapse model as 



We use the following expressions (Kitayama & Suto 1996) for the critical over-density A n i: 



Function 



Our present method to compute the PVDF is schematically shown in Figure [TJ. We 
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Fig. 1. Schematic illustration of our procedure to compute the pairwise velocity distribution function. 
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where uj vir and 77 vir are respectively given as uj y - 1t = l/fi V ir — 1 an d r) vil = cosh~ 1 (2/f2 vir — 1), in 
terms of the density parameter at the collapse time, Q vir . 

In practice, we focus on three specific profiles, (i) the original NFW profile with a = 1 
and v = 2 (Navarro et al. 1997), (ii) the modified NFW profile with a = u = 3/2 indicated by 
higher-resolution simulations (Fukushige & Makino 1997, 2001a,b; Moore et al. 1998; Jing & 
Suto 2000), (iii) the Hernquist profile with a = 1 and v = 3 (Hernquist 1990) for which the 
analytic expression of the phase-space distribution function is known. 

The two parameters r s and r v ; r are not independent, and are related in terms of the 
concentration parameter: 

r viT (M,z) 



c(M,z) 



(7) 



r s {M,z) " 

The condition that the total mass inside r vlv is equal to M relates 5 C to c. Therefore the 
halo mass-dependence of the above profiles is entirely specified by c = c(M). In the case of 
the original NFW profile, we use the approximate fitting function from the simulation data of 
Bullock et al. (2001): 

8.0 f M \"°' 13 

CB (M) = — (^t^H • (8) 
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For the other profiles, we first compute the amplitude of the two-point correlation functions 
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of dark matter following the procedure of Seljak (2000) and Ma & Fry (2000), and then find 
the amplitude of the concentration parameter which reproduces the Peacock - Dodds (1996) 
fitting formula. This calibration yield c(M) = cb(M)/2 for the modified NFW profile (Oguri 
et al. 2001), and c(M) = cb(M)/3 for the Hernquist profile, which we adopt throughout the 
analysis below. 

2.2. Phase-space distribution function in a halo 

Our next task is to compute the phase-space distribution function in a single halo from 
the given density profile (f|). While Sheth (1996) and Sheth et al. (2001) simply adopt the 
Gaussian velocity distribution function, we eliminate this assumption and derive the velocity 
distribution function in a fully consistent manner. 

For this purpose, we make use of the Jeans theorem which states that for a spherically 
symmetric and stationary system, the solution of the collisionless Boltzmann equation can be 
expressed as a function of the specific binding energy, E = ip(r) —v 2 /2, alone. Here we define ip 
as the minus of the gravitational potential satisfying the boundary condition of tp(r — > oo) — > 0. 

One may wonder whether halos in hierarchical universes that should experience re- 
peated merger and destruction continually are well approximated as stationary. Nevertheless 
Natarajan, Hjorth & van Kampen (1997) and Hanyu & Habe (2001) found that the phase-space 
distribution function directly estimated from their particle simulations agrees well with that 
derived from the Jeans theorem. Thus the above assumption is justified, at least empirically. 

Then the phase-space distribution function F(E\M) in a single halo is directly computed 
from its given density profile p(r\M) as follows (e.g., Binney & Tremaine 1987): 



Figure ^| plots the dimensionless phase-space distribution function f(e) = 
F(E\M)(Gm s /r s ) 3 / 2 /(5 c p) evaluated numerically from equation (^), where m s = Anr^5 c p 
is the characteristic mass of the halo and e = Er s j \Gm s ) is the dimensionless specific energy. 
The Hernquist model has an analytical solution for F(E\M) which is reproduced by our 
numerical result almost within an accuracy of 2% except for the tiny region e ~ 0, where the 
error reaches at 7% but the effect is safely negligible for later analysis. Since the three halo 
profiles that we adopt have a central cusp, f(e) diverges at a corresponding value of e. The 
modified NFW profile has f(e) which extends more broadly up to e ~ 2 reflecting the stronger 
central concentration than that of the original NFW case. 

2. 3. Single-particle velocity distribution function in a halo 

Once the phase-space distribution function is given, one can also compute one- 
dimensional single-particle velocity distribution function along a particular direction by in- 
tegrating over the other two components. Assuming the isotropic velocity distribution, one 
has 




(9) 
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Fig. 2. The dimensionless phase-space distribution function /(e) = F(E\M)(Gm s /r s ) 3 ^ 2 /(S c p) as a func- 
tion of dimensionless binding energy e = Er s / (Gm s ); NFW(q = 1) profile (solid), NFW(a = 1.5) profile 
(dotted) and the Hernquist profile (dashed). 
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where we project along the direction of V\ and the quantity E\ = ip — vf/2 is the corresponding 
binding energy. Figure |3| shows the dimensionless velocity distribution function v s fi(v\M;r) at 
r/r s = 1 (thin lines) and r/r s = 10 (thick lines), where v s = (Gm s /r s ) l l 2 is the scaling velocity. 
The figure indicates that the one-dimensional velocity distribution function in a halo can be 
reasonably approximated by the Gaussian: 

f 1 (v\M:r)= ,— 1 expf ^ — - ) , (11) 

J K ' ' J V2^a(r\M) y \ 2a 2 (r\M) J V ; 

although it has a sharp cutoff around the escape velocity of the halo, v csc = (2ip) 1 / 2 . 

In Figure |3], the dashed lines indicate the Gaussian-fit which has the same velocity 
dispersion evaluated from equation ([TP]) . It seems that the empirical Gaussian approximation 
can reasonably reproduce the PVDF, and thus we use the approximation in the numerical 
integrations below so as to reduce the computational time. 

Figure |] plots the velocity dispersion a(r\M) computed from the best-fit Gaussian, 
which clearly shows the scale-dependence that was neglected in the previous analysis (Sheth 
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Fig. 3. Single-particle velocity distribution function in a single halo at r/r s = 1 (thin-solid) and r/r a = 10 
(thick-solid). The dashed lines show the corresponding Gaussian fits (see eq.JTl]])- Left: NFW(a = i); 
Middle: NFW(a = 1.5); Right: Hernquist profile. 

1996; Sheth et al. 2001). Note also the different scale-dependence from the circular velocity 
V c = Gm(r)/r, where m(r) is the mass inside the radius r (thin lines in Fig.|j). In the subsequent 
modeling of the PVDF, we use the Gaussian approximation with the fitted o~(r\M) rather than 
repeating the full numerical integration. 

2-4- Pairwise relative velocity distribution function 

Finally we are in a position to estimate the PVDF combining the above results. Since 
we are interested in small scales, the particle pairs with the corresponding separations are 
approximated to reside in the common halo. Then Sheth (1996) derived the following expression 
for the PVDF: 

f(vi2,r 12 ) =-Af -1 / dMn(M) J d 3 ri d 3 r 2 p( ri \M)p(r 2 \M) 

x J dv 1 dv 2 fi(v 1 \M;r 1 )f 1 (v2\M]r2)S D (r 12 -\r 1 -r2\)S D (v 12 -v 1 +v 2 ), (12) 

where r 12 is the pair-separation and M is the normalization factor given by 

M = jdMn(M) |c/ 3 r 1 rf 3 r 2 p(r 1 |M)p(r 2 |M)5 D (r 12 -|r 1 -r 2 |). (13) 

The relation between the quantities used in the expressions ( |I"2"D and ([T3| ) are schematically 
summarized in Figure ||. 

While Sheth (1996) adopted the singular isothermal sphere p(r) oc r~ 2 and therefore 
the Maxwellian for j\(v\M]r) with r-independent velocity dispersion, we are able to evaluate 
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Fig. 4. The peculiar velocity dispersion a(r\M) of dark matter particles in a single halo resulting from 
the Gaussian fit as a function of position r/r s (thick lines). For comparison, the thin lines show the 
circular velocity V c evaluated from the relation V c (r) = Gm(r)/r, where m(r) represents the mass inside 
the radius r. Solid line, dotted line, and dashed lines represent the results in NFW(a = 1), NFW(a = 1.5), 
and Hernquist profiles, respectively 




Observer 



Fig. 5. The schematic picture of the relations of the relative positions and velocities used in the expressions 
(12) and (|l3|). The large circle and the two small circles denote a halo and a pair of particles, respectively. 
The separation length of the pair r 12 is defined by \ri — r 2 | and the relative pairwise velocity w 12 is defined 
by vi — V2- The quantities v\ and V2 are the line-of-sight component of the three dimensional velocities iti 
and i6 2 , respectively. 
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equation fll2|) in more realistic situations as described in the preceding subsections. 

Applying the Gaussian fit (eq. fill) for fx(v\M;r), equation fljjQ is rewritten as 

h,(v„\r»)=N- 1 fdM„(M) ItrJfr* ^\M) ^\M) 



x J dv\dv2 exp 
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2ira(r 1 \M)a(r 2 \M) 

2a*( ri \M) 2a\r 2 \M)\ ^ ~ ^ ~ ^ ^ ~ Vl + %h){U) 
The integrals over the two velocity components V\ and v 2 , and also over the directions of the 
position vectors can be performed analytically, and equation (|14|) reduces to 

fi2(v 12 \ri2) = 2 _ 2 / dMn(M)r 12 dr 2 dr x 

^ r \2p " ,M min(^l2) Jmax(0,ri2 -r s c) ^Ir^-r^ 



p(n\M) p{r 2 \M) 
x rir 2 exp 

^2n{a 2 ( ri \M)+a 2 (r 2 \M)} 
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2{a 2 (r 1 |M) + <x 2 (r 2 |M)} 



(15) 



where M min (r 12 ) is the minimum mass of the halo including the pair with separation r 12 (i.e., 
r vir (M min ) > r 12 /2), and the normalization factor M is now given by 

\ roo rr B c rmin(r B c,ri2+r2) 

M=- — — - / dMn(M)r 12 dr 2 dnnr^n^p^M). (16) 

4VTr 12 p "'M m in(fi2) Jmax(0,ri2 — r s c) J|ri 2 -f2 

Actually it turns out that this term corresponds to the one-halo contribution of the two-point 
correlation function, £ih(fi2), m the dark halo approach (Seljak 2000, Ma & Fry 2000). 

Figure ^| plots the resulting PVDF (in the LCDM model) at r 12 = lh~ l , 0.3/i _1 and 
0.1/i~ 1 Mpc against the pairwise velocity v\ 2 normalized by the PVD at each separation. We 
adopt the Press-Schechter mass function for definiteness. Note that the quantitatively similar 
behavior is obtained for other cosmological models, but with different PVDs. The PVDF 
for small separation pairs ir\ 2 <^ l/i _1 Mpc) is well described by the exponential distribution. 
As r\ 2 increases, the central region resembles the Gaussian distribution while the exponential 
tail is still clear at large velocities. Neither the inner nor outer slope of the density profile 
produces any systematic difference in the non-Gaussian tails of PVDF, in contrast to the single- 
particle velocity distribution in Figure |3|. This qualitative behavior is in complete agreement 
with the result of Sheth (1996) assuming the scale-free model and the singular isothermal 
sphere. Therefore we conclude that the exponential distribution of the PVDF is a rather 
general consequence in the gravitational instability picture fairly independent of the underlying 
cosmological model. 

3. Pairwise relative peculiar velocity dispersions 

3.1. Dark matter particles 

We have shown that the shape of the PVDF is well approximated by the exponential in a 
fairly insensitive manner to either the cosmological model or the dark halo density profile. Note, 
however, that this does not implies the PVD is determined independently of the cosmology, but 
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Fig. 6. Pairwisc velocity distribution function averaged over the halo mass function in the LCDM uni- 
verse. Left: NFW(a = 1); Middle: NFW(a = 1.5); Right: Hernquist profile. Solid, dotted and dashed lines 
indicate the results at the pair separation of l/i -1 Mpc, 0.3/i _1 Mpc, and 0.1/i _1 Mpc. 

rather it is the basic source for the cosmological model-dependence, through the mass function 
of halos. Since the PVDF is already obtained, one may evaluate the PVD o"i 2 (ri 2 ) directly as 



^l 2 0l2) 



dv 12 fi 2 (v 12 ;r 12 ) v 



12- 



(17) 



We note, however, that the above expression is valid only when the one-halo term £ih(Vi2) is 
sufficiently larger than unity. If one takes account of particle pairs residing in two different halos 
that we neglect in the present modeling, one should rather replace the normalization factor N 
by 1 + £(7*12) since the factor physically corresponds to the relative probability of finding a pair 
at separation r 12 . This consideration implies the normalization of the PVD from the one-halo 
contribution should be 



'12 



(ria) 



£ih(ri 2 ) f°° 



dv 12 fi 2 (v i 2 ;ri 2 ) v\ 2 . 



(18) 



l + £(ri 2 ) 

In fact, this agrees with equation (21) of Sheth et al. (2001). We compute the two-point 
correlation function, £(7*12) = £ih(^i2) +£ 2 h(^i2), on the basis of the dark halo approach (Seljak 
2000), and we also confirm that the resulting £(r 12 ) is in good agreement with the fitting formula 
of Peacock & Dodds (1996) by an appropriate choice of the concentration parameter c(M) as 
discussed in subsection [2~T| . 

Figure shows the PVD calculated according to equation (|TJ). The comparison among 



10 



7] 



1200 
1000 
800 
600 



~ 400 



200 



1 Mill 1 III 

NFW(«=1) 

/A 


i i i i i in i ii 

NFW(a = 1.5) 


1 1 1 1 1 1 1 1 III 

Hernquist 




\v\ - 




X \ * 

\\ - 


_^ . ,i _ 




SCDM \ \ \, 


\\\- 




LCDM \ \ \ 






OCDM \ \ \ 

- LCDM with \ \ \ - 

Maxwellian \ \ \ 


\\\- 

1 I No. -4- 


i i iM-k 





0.1 1 0.1 1 0.1 1 

r 12 [h _1 Mpc] r 12 [h -1 Mpc] r 12 [h _1 Mpc] 

Fig. 7. The pairwisc velocity dispersion of dark matter particles. Left: NFW(q = 1); Middle: 
NFW(a = 1.5); Right: Hernquist profile. Solid, dotted and dashed lines indicate the results in SCDM, 
LCDM, and OCDM. Long-dashed line shows the result in LCDM under the assumption that each halo 
has the Maxwellian velocity distribution function with constant velocity dispersion. 

the different model predictions indicates that the amplitude of PVD sensitively depends on 
the cosmological parameters through the mass function n(M), but is almost insensitive to the 
density profile of dark halo. Note, however, that this is partly because we have chosen the value 
of c(M) so as to reproduce the same £(r) irrespectively of the density profile. We also show 
the result with neglecting the scale-dependence of the velocity distribution in the case of the 
LCDM model (long-dashed lines). This indicates that the isothermal approximation (Sheth et 
al. 2001) is quite acceptable in predicting the PVD. 

We attempted an empirical fitting to our numerical results at different redshifts (Fig.|8|) 
by adopting a power-law form: 

t> 



012 (Z12) = A 



X12 



(19) 



in terms of the comoving pair separation x 12 = r 12 (1 + z) . The values of the amplitude A and the 
power-law index p fitted for 0.01/z _1 Mpc < X12 < i m ax are listed in Table |l] in the LCDM model, 
where x max is the comoving separation at which the PVD becomes maximum. Since the two- 
halo term ^2h(^i2) becomes dominant contribution to £(Vi 2 ) for larger separations, equation ([TP] ) 
becomes inaccurate for x > x max . The fit is accurate within 5 %, which is comparable to the other 
systematic errors including the numerical integration or the Gaussian approximation. This 



11 



x 



1200 

1000 

800 

600 

~~ 400 
b 

200 



0.01 0.1 1 0.01 0.1 1 0.01 0.1 1 
x 12 [h _1 Mpc] x 12 [h _1 Mpc] x 12 [h _1 Mpc] 

Fig. 8. The pairwise velocity dispersion of dark matter particles in LCDM at z = 0, 0.5, 1.0, 1.5, 2.0 and 
2.5 from top to bottom. Left: NFW(a = 1); Middle: NFW(a = 1.5); Right: Hernquist profile. 

result may be compared with the independent prediction based on the cosmic virial theorem 
(Peebles 1976; Suto 1993; Suto & Jing 1997). For reference, if the correlation function of dark 
matter is given by (r/5.4/i _1 Mpc) -1 ' 8 , the cosmic virial theorem implies that 
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(20) 



\ v l/i~ 1 Mpc / 

where Q is the normalized amplitude of the three-point correlation function. The value of Q 
for dark matter clustering is somewhat uncertain, but may be close to 2 in LCDM (Suto 1993), 
and thus this estimate is fairly consistent with our prediction presented here. 

Table 1. Power-law fits to the relative peculiar velocity dispersion in the LCDM model. 
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3.2. Effect of biasing 

It is well known that the PVD of the observed galaxies is generally smaller than the 
value predicted in current popular models (e.g., Davis & Peebles 1983; Mo, Jing & Borner 
1993; Suto 1993; Suto & Jing 1997). This may be interpreted as a manifestation of the spatial 
biasing of galaxies relative to dark matter. 

Jing, Mo & Borner (1998) analyzed the Las Campanas Redshift Survey galaxies and de- 
veloped a phenomenological biasing model, CLuster underWeight bias (CLW, hereafter) which 
successfully accounts for the amplitudes of the two-point correlation function and the PVD 
simultaneously. More recently Jing, Borner & Suto (2002) performed the similar analysis of 
galaxies in the PSCz catalog (Saunders et al. 2000) which are selected from the InfraRed 
Astronomical Satellite (IRAS) Point Source Catalog (PSC; Beichman et al. 1988). They found 
that the IRAS selected galaxies, which are likely to be dominated by late-types, have signifi- 
cantly smaller PVD than those in other catalogues. In addition, they applied the CLW bias 
scheme to mock samples from iV-body simulations and concluded that the PVD of the PSCz 
galaxies is significantly smaller than those predicted from the CLW bias in the popular CDM 
models. In this subsection, we revisit this issue combing the biasing effect with our analytical 
model of the PVD in a complementary manner to the direct method using the iV-body data. 

The CLW scheme can be applied to our analytical model easily. According to Jing et 
al. (2002), we adopt that the selection function of galaxies in a halo of mass M is proportional 
to (M/1O 14 M ) _/3 . This is simply equivalent to replacing the mass function n(M) in equation 
(|15"[) by n(M)(M/lO 14 M ) _/3 . Because this biasing model puts lower weight on the massive 
clusters where the velocity dispersion is large, increasing (3 suppresses the mean PVD. In fact, 
the Las Campanas redshift survey data are consistent with (3 = 0.08, and the PSCz data prefer a 
much larger value (3 = 0.25. Physically speaking, this phenomenological prescription should be 
understood as the dependence of the efficiency of galaxy formation on the mass of the hosting 
halo. 

In addition we consider another biasing to take into account the observed density- 
morphology relation of galaxies. Since spiral galaxies preferentially avoid the central region of 
massive clusters (i.e., halos in the present context), the PVD of spirals is generally suppressed 
especially in the case of the modified NFW profile that has stronger central concentration (c.f., 
Fig.|]). We attempt to incorporate this effect by introducing the selection probability that 
depends on a distance from the center of the halo: 

^ M) - a+ dw) k (21) 

We set a = 0.2 and b = 0.6 so as to reproduce the observed ratio of spirals to ellipticals; 2:8 in 
the inner part and 8:2 in the outer part. In practice, we calculate the PVD by replacing p(r\M) 
by p(r\M)p(r\M) in equation fllB]) . 

The results are shown in Figure || for the modified NFW profile in the LCDM model. 
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Fig. 9. Pairwisc peculiar velocity dispersion for galaxies with empirical biasing scheme in LCDM; Solid: 
dark matter particles. Dotted: CLW with j3 = 0.08, Short-dashed: CLW with j3 = 0.25, Long-dashed: CLW 
with (3 = 0.25 and density-morphology relation. The filled squares indicate the values for the PSCz galaxies 
estimated by Jing et al. (2002). 

The dotted and short-dashed lines represent the results taking account of the CLW bias effect 
while the long-dashed lines consider the density-morphology relation (eq. ||21||) in addition to 
the CLW bias. The degree of suppression of the PVD is in agreement with the simulation 
results of Jing et al. (2002). Even with the density-morphology relation, the PVD of the IRAS 
PSCz galaxies is too small to be reconciled in the current model as Jing et al. (2002) claimed. 

4. Summary 

We have presented a detailed prediction for the pairwise peculiar velocity distribution 
function (PVDF) applying the dark matter halo approach. In particular, we have derived the 
PVDF in a direct and self-consistent manner with the assumed density profile for dark matter 
halo for the first time. On the other hand, we neglect the halo-halo contribution assuming that 
any pair of particles resides in a common halo. Thus our predictions are quantitatively valid 
only on small scales < l/i _1 Mpc, but our result turns out to be fairly close to the previous one 
by Sheth (1996) and Sheth et al. (2001) who assumed the isothermal velocity distribution in a 
single halo. In this sense, our independent approach may be regarded as to provide an empirical 
justification for their simplifying assumptions, and also our predictions are fairly accurate on 
those small scales. 

We have shown that the shape of the PVDF is well approximated by the exponential in 
a fairly insensitive manner to either the cosmological model or the dark halo density profile. 
The dependence on the PVD on the halo density profiles that we employed is also fairly small, 
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yielding the difference less than about 10 percent. 

We have also obtained a practical fitting formula for the PVD of dark matter particles 
(eq. ||19|| ) at different cosmological models as a function of the pair separation. This may be 
useful in modeling redshift-space distortion of clustering. The result is in reasonable agreement 
with the estimate on the basis of the cosmic virial theorem. Furthermore we apply an empirical 
biasing scheme into our model and attempt to predict the PVD of galaxies. We can reproduce 
the previous simulation results on the basis of our analytical method, and also confirmed that 
the very small PVD estimated for the PSCz galaxies ( Jing et al. 2002) is difficult to be reconciled 
with a simplistic biasing model and/or the underlying CDM model. 

The discrepancy between the prediction and the observation shown in section |3.2| in- 
dicates the presence of velocity bias, in addition to the other selection effects. In fact, each 
luminous galaxy is a clump composed of the baryon as well as dark matter particles, whose bulk 
motion might not trace the random motion of the individual dark matter particles. In this case, 
the galaxy-galaxy interaction through tidal field or gas pressure of baryon might be an impor- 
tant source for velocity bias. At least, our present treatment using dark matter halo approach 
can provide a quantitative prediction for the PVD of the dark matter particles. Therefore next 
step, the effect of velocity bias including these interactions should be incorporated into our 
scheme to account for the PVD of the galaxies. 
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